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Abstract We report a new precise determination of the leptonic widths of the J/psi meson performed with 
the KEDR detector at the VEPP-4M e + e~ collider. The measured values of the J/psi parameters are: 

r ee xr ee /r =0.3323±0.0064(stat.) ± 0.0048 (syst.) keV, 
R e xr MM /r = 0.3318±0.0052(stat.) ±0.0063 (syst.) keV. 

Assuming e/i universality and using the table value of the branching ratios the leptonic F« = 5.59±0.12keV 
width and the total F = 94.1 ± 2.7keV widths were obtained. We also discuss in detail a method to calculate 
radiative corrections at a narrow resonance. 
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1 Introduction 



The J/tp meson is frequently referred to as a hy- 
drogen atom for QCD. The electron widths T ce of 
charmonium states are rather well predicted by po- 
tential models [l|, Q . The uncertainty in the QCD lat- 
tice calculations of T ee gradually approaches the ex- 
perimental errors Q. The full and dileptonic widths 
of a hadronic resonance, T and T«, describe funda- 
mental properties of the strong potential 

In this report we discuss the results of the J/ip 
meson observation in leptonic decay channels. Study 
of the e + e~ — > J/ip — > £ + £~ cross section as function 
of energy allows one to determine the leptonic width 
T u and its product to the decay ratio T ee xT u /T thus 
the total width T can be also found. The product 
r ee x T u /T determines the peak cross section while 
the leptonic width T u is contained in the interference 
wave magnitude. Due to smallness of the interference 
effect the experimental accuracy of the Y u determina- 
tion is rather poor. However, the branching ratio B u 
is known with the accuracy of 0.7% from the cascade 
decay ip(2S) — > J/ipir + Tr~ thus we report the high 
precision results on T ee x T ee /T and T ee x T^/T and 
use the T u value to check the analysis consistency 
only. 



The extraction of resonance parameters from the 
measured cross section requires the accurate account- 
ing of radiative corrections. The Sec. 8.2.4 of the 
highly cited report Q treats the radiative corrections 
to e + e~ — > J lij}— > £ + £~ cross section in the way con- 
tradicting to that used in the experiments fig and 
our work 0] therefore me start with the discussion of 
this issue. 

2 Radiative corrections to J/ip pro- 
duction and decays 

In virtually all experimental analyses it is assumed 
that the resonant contribution to the cross section of 
e + e~ — > J/ijj — > £ + £~ is proportional to the product 
r ee x T u I r where r ee and Y u are so called experi- 
mental partial widths [8j recommended to use by the 
Particle Data Group since 1990: 



r« = Bu(n~t) x r : 



i-n 



(i) 



where Bu^ ni ) is the branching ratio as it is measuring 
experimentally, T° e is the lowest order QED partial 
width and n is the vacuum polarization operator ex- 
cluding J/ip contribution. In contrast, the Sec. 8.2.4 
of Ref. [4[ proposes that the resonant contribution is 
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proportional to T ee x T ( ^_/T = T° ee x T u /T. 

According to Ref. [9] the cross section of the 
single-photon annihilation e + e~ — >£ + £~ can be writ- 
ten in the form 



= / dx ■ 



a ((l-x)s) 
l-n((l-a:)s) 



■f(s,x), 



(2) 



where the f(s,x) is calculated with a high accuracy, 
the II(s) represents the vacuum polarization opera- 
tor and tr (s) in the Born level cross section of the 
process. 

Assuming the Breit-Wigner shape for <j 



i27iT° r° 



{s-M 2 ) 2 + M 2 T 2 



(3) 



and replacing II(s) with n mentioned above, one re- 
produces the result of the Sec. 8.2.4 of Ref. [4|. 

However, the Born level cross section of the 
e + e~ — > £ + £~ process is the smooth function of s 
therefore the resonance behavior of the cross section 
^) is due to energy dependence of the full vacuum 
polarization operator II containing the resonant con- 
tributiorQ- One has 13 = n + TIr with nonresonant 

n = n ee +n MAl +n TT +n 9? and 



n fl (s) = 



3r f ' s 



i 



a M o s-M$+iM o r 



(4) 



where M , T and T^J are the "bare" resonance mass 
and widths. 

The formula ((2]) gives the cross section without 
separation to the continuum, resonant and interfer- 
ence parts. To obtain the contribution of the reso- 
nance, the continuum one must be subtracted from 
the amplitude. It can be done with the equality 



1 



i-n -n fl ( s ) 
i i 



3r° 



l-IIo (l-IIo) 2 a M 0S -M 2 + iMT 
in which both M and T depend on s: 



(5) 



M' 2 = Ml + 



Mr = M o r - 



^^Re^ 
a M l-IIo ' 

a M 1 



n 



(G) 



In a vicinity of a narrow resonance this dependence is 
negligible thus the resonant contribution can be de- 
scribed with the Breit-Wigner amplitude containing 
"dressed" parameters M«M(M 2 ), r«f(M 2 ). Due 
to the extra power of the vacuum polarization factor 
1/ 1 1 — n 1 m the second term of ([5]) the resonant part 



of the e + e" — > £ + £~ cross section is proportional to 
r ee x r«/r and does not depend on T^J explicitly. 

The analytical expressions for the e + e~ — > £ + £~ 
cross section in the soft photon approximation were 
first derived by Ya. A. Azimov et al. in 1975 10 1. 
With some up-today modifications one obtains in the 
vicinity of a narrow resonance 

'da} ee " MM _ f da } ee " MM _ 3 



dfl 



TJH +— -(l + 4 f )(l + cos^)x 



QED 



4M 



3rJ ""-ImJ- 2a ^ 



r 



YM 



M 



Re- 



T 



l-IIo 



(7) 

where a correction S S f follows from the structure func- 
tion approach of 



and 



with 



7T/3 



37 



1 



W 



96 12 36 



m, 



M/2 



sinTr^ \-W + M-iY/2 



7T \ Tflp 2 



(8) 
(9) 
(10) 



The terms proportional to ImJ 7 and Re J 7 describe 
the contribution of the resonance and the interference 
effect, respectively. 

Originally in Ref. [l(J the electron loops only were 
taken into account in n while the terms < /3 2 were 
omitted including the 7r/3/sin7r/3 factor ll| in 
For the e + e~ final state one has 
daY e " ee ^ (da" 

dn) ~[dn. 



i ( ■< c 



+ 



QED 



i f 9 r 2 
— -— (i 

M 2 I 4 YM 



-cos 2 9) il- 



ia r ee 
T m" 



(i- 



-COS V) — 



(1 



-5 st ) imJ 7 - 
■f cost?) 2 



(l-cos0) 



Re J 7 



(11) 

where the relative accuracy of the interference term 
is about (3 (7.6% for J/ip). That is sufficient for the 
analysis reported. 

For the nonresonant contribution <7q B d the calcu- 



lations of [12|, [13| can be used implemented in the 



event generators BHWIDE [14j and MCGPJ [15 1. 

In order to compare the theoretical cross sec- 
tions ([7|) and (fTT|) with experimental data, it is neces- 
sary to perform their convolution with a distribution 
of the total collision energy which is assumed to be 
Gaussian with an energy spread a w : 

1 ( (W-W ) 2 * 



p(W) = 



2tt<j w 



exp 



*We are grateful to V. S. Fadin for clarification of this issue. 
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where W is an average cm. collision energy. 

3 VEPP-4M collider and KEDR de- 
tector 



cesses. The on-line luminosity is measured by two 
independent single bremsstrahlung monitors. 

4 Experiment description 



The VEPP-4M collider [16| can operate in the 
wide range of beam energy from 1 to 6 GeV. The 
peak luminosity in the J/ip energy region is about 2x 
lO^cm-V 1 . 




CO, laser, A=10.59 fj.m 

Fig. 1. VEPP-4M/KEDR complex with the 
resonant depolarization and the infrared light 
Compton backscattering facilities. 

One of the main features of the VEPP-4M is a 
possibility of precise energy determination. The reso- 
nant depolarization method 17, 3] was implemented 



at VEPP-4 from the beginning of experiments in early 
eighties for the measurements of the J/if) and tp(2S) 
mass with the OLYA detector and T family mass 



with the MD-1 [19j detector. 

At VEPP-4M the accuracy of the energy calibra- 
tion with the resonant depolarization is improved to 
about 10~ 6 . The interpolation of energy between cal- 
ibrations in the J/tp region has the accuracy of 
6-1CT 6 (-10 keV). 

In 2005 a new technique developed at the BESSY- 
I and BESSY-II synchrotron radiation sources 2 ll 1 2 2 1 



was adopted for VEPP-4M. It employs the infrared 
light Compton backscattering and has a worse preci- 
sion (50-^70 keV in the J/ip region) but, unlike the 
resonant depolarization, can be used during data tak- 
ing. 

The KEDR detector [23| includes the vertex de- 
tector, the drift chamber, the scintillation time-of- 
flight counters, the aerogel Cherenkov counters, the 
barrel liquid krypton calorimeter, the endcap Csl 
calorimeter, and the muon system built in the yoke 
of a superconducting coil generating a field of 0.65 
T. The detector also includes the scattered electron 
tagging system for studying of the two-photon pro- 



In April 2005, the 11-point scan of the J/tp) has 
been performed with the integral luminosity of 230 
nb _1 . This corresponds approximately to 15000 
J/tp — » e + e" decays. During this time, 26 calibrations 
of the beam energy were done using the resonance- 
depolarization method. 

Single bremsstrahlung and Bhabha scattering to 
the endcap calorimeter were used in the relative mea- 
surement of luminosity. The absolute calibration of 
the luminosity was performed using the large angle 
Bhabha scattering in the F ee x F ee /r analysis. 

fobs (nb) 
1600F 




3088 



3092 



3096 



3100 



3104 3108 

E w (MeV) 



Fig. 2. Observed e + e — ► hadrons cross section 
according to the results of the J/tp scan. 

Figure [2] shows the observed e + e~ — > hadrons cross 
section at the J/tp energy region. These data were 
used to fix the resonance peak position and to de- 
termine the beam energy spread. The value of the 
J /ip mass agrees with the earlier VEPP-4M/KEDR 
experiments [io| . 

5 Data analysis 

In our analysis we employed the simplest selection 
criteria that ensured a sufficient suppression of multi- 
hadron events and the cosmic-ray background, please 
see Ref. 



24] for details. 



In order to measure the resonance parameters in 
the e + e~ channel, the set of events was divided into 
ten equal angular intervals from 40° to 140°. At the 
i-th point in energy Ei and the j-th angular interval 
9j , the expected number of events was parameterized 
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as 



N mqf {E i ,9 j )=H c y.C{E i )x 



a 



a 



EShaMiaV®*' ®j) 



habha 



(12) 



where C(Ei) is the integrated luminosity measured 
by luminosity monitor at the i-th point; crJg S eor , cr*^ t e e ° r 
and tJg 1 5 1 e ° 1 ^ ha are the theoretical cross sections respec- 
tively for resonance, interference and Bhabha contri- 
butions; ef™, e?X and e^ abha are detector efficien- 
cies obtained from simulated data. 

In this formula the following free parameters were 
used: 

1. the product r ec x r ec /r, which determines the 
magnitude of the resonance signal; 

2. the electron width r ee , which specifies the am- 
plitude of the interference wave; 

3. the coefficient 7Z C , which provides the absolute 
calibration of the luminosity monitor. 





* . A, 


10° < 6 


< 65° 
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MeV 
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65° SS < 90' 
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Fig. 3. Fits to experimental data for e + e~ — > 
e + e~ process at J/tp energy region for four 
angular ranges. 

We note that the coefficient IZc partially accounts 
the possible difference between the actual detection 
efficiency and simulation in the case where these dif- 
ference do not depend on the scattering angle or the 
beam energy (or the data taking time) thus the sub- 
stantial cancellation of errors occurs. 

Figure [3] shows our fit to the data for four angular 
intervals. The joined fit in ten equal intervals from 
40° to 140° produce the following basic result: 

T ee x r ee / T = 0.3323 ± 0.0064 (stat.)keV, 
7e £ = 93.4±0.7(stat.)%, 

r ee = 5.7±0.6(stat.)keV. (13) 

Due to different angular distributions for Bhabha 
scattering and resonance events, subdivision of the 



data into several angular bins decreases a statistical 
error for T ec xT ee /T by 40-=- 50%. The electron width 
obtained by the fit has a statistical error of about 10% 
and agrees with the world-average value. 




3088 3090 3092 3094 3096 3098 3100 3102 3104 3106 3108 

W, MeV 

Fig. 4. Fit to experimental data for e + e~ — > 
H + fiT process at J/tf) energy region. 

Similarly to (TT2")) . the expected number of e + e~ 
\i~ events was parameterized in the form: 



AT exp (£i)=^x£(£i)x 



'bg 



<(E i )-e«Zr(E i )+ 



(14) 



. y T 

ale 



with the same meaning of TZc an d C(Ei) as in (TT! 
C{Ei) is multiplied by the sum of the products of 
theoretical cross sections for resonance, interference 
and QED background and detection efficiencies as 
obtained from simulated data. IZc w &s fixed by re- 
sult (TT5]) . T t is the live data taking time. Unlike (fT2"j) 
there is only one angular interval from 40° to 140°. 
The following free parameters were used: 

1. the product T ee x r MJJ /r, which determines the 
magnitude of the resonance signal; 

2. the square root of electron and muon widths 
y/T ee T MM , which specifies the amplitude of the 
interference wave; 



3. the cosmic events rate F c( 
tion criteria for the e + e~ - 



nic passed the selec- 
■ H~ events. 

Due to variations of luminosity during the experiment 
it is possible to separate cosmic events contribution 
(-Pcosmic ' Ti) from nonresonant background contribu- 
tion « cor (^)-eb g m (^)-^)- 

Figure |4] shows our fit to the e + e" — > fi + \i~ data. 
It yields the following result: 



T ee x r MM /r = 0.3318±0.0052 (stat.)keV, 



x/T ee x T MM = 5.6 ± 0.7 (stat.) keV. 



(15) 
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As can be seen from (|15j) the statistical error of 

r ee xr MM /r is about 1.6%. 

6 Discussion of systematic uncer- 
tainty 

The most significant systematic uncertainties in 
the r ee x r ee /r and T ee x T^/T measurements are 
listed in Tables Q] and [2j respectively. 



Table 1. Systematic uncertainties in r ee 


xr ee /r. 


Systematic uncertainty source 




Error % 


Luminosity monitor instability 




0.8 


Offline event selection 




0.7 


Trigger efficiency 




0.5 


Energy spread accuracy 




0.2 


Beam energy measurement (10-30 kcV) 




0.3 


Fiducial volume cut 




0.2 


Calculation of radiative correction 




0.2 


Cross section for Bhabha (MC generators) 




0.4 


Uncertainty in the final state radiation (PHOTOS) 


0.4 


Background from J/ip decays 




0.2 


Fitting procedure 




0.2 


Quadratic sum 




1.4 


Table 2. Systematic uncertainties in r ee 


xr w /r. 


Systematic uncertainty source 




Error % 


Luminosity monitor instability 




0.8 


Absolute luminosity calibration by e + e — data 




1.2 


Trigger efficiency 




0.5 


Energy spread accuracy 




0.4 


Beam energy measurement (10-30 kcV) 




0.5 


Fiducial volume cut 




0.2 


Calculation of radiative correction 




0.1 


Uncertainty in the final state radiation (PHOTOS) 


0.5 


Nonresonant background 




0.1 


Background from J/ip decays 




0.6 


Quadratic sum 




1.8 



A rather large uncertainty of 0.8% common for the 
electron and muon channels is due to the luminosity 
monitor instability. It was estimated from compar- 
ing the results obtained using the on-line luminosity 
of the single bremsstrahlung monitor and the off-line 
luminosity measured by the e + e~ scattering in the 
endcap calorimeter. 

The essential source of uncertainty is an imper- 
fection of the detector response simulation resulting 
in the errors in the trigger and offline event selection 
efficiencies. It was studied using collected data and 
the correction of 0.75 ±0.7% was applied. 



The dominant uncertainty of the T ee x T^/T re- 
sult is associated with the absolute luminosity cal- 
ibration done in e + e~-channel. It includes the ac- 
curacy of the Bhabha event generators, the statistic 
error of IZc parameter (|13[) and the residual (after 
correction using simulated data) efficiency difference 
for e + e~ and pi + ji~ events. The additional correction 
applied to this difference is — 0.5 ±0.9%. 

The other sources of uncertainty are discussed in 
Ref. [24 1 . 



7 Results and Conclusion 

The new measurement of the T ee x T ee /T and 
r ee xr MM /r has been performed at the VEPP-4M col- 
lider using the KEDR detector. The following results 
have been obtained (in keV): 

r ee xr ee /r =0.3323±0.0064(stat.)±0.0048(syst.) 
r ee xr MM /r = 0.3318±0.0052(stat.) ±0.0063 (syst.) 



Tee X ^fjp/ r 



FRAM 1975 ■ ■ — 

DASP 1975 

BaBar 2004 — 

CLEO-c 2006 !• 

KEDR 2009 ~ 

Tee X Tee/ I 



SPEC 1975 
FRAM 1975 
FRAG 1975 
DASP 1979 
KEDR 2009 



0.2 



0.3 



0.4 



0.5 



~0J) keV 



Fig. 5. Comparison of r ee X r ee /r and r ee x 
r MM /r measured at different experiments 
mentioned in with KEDR 2009 results. 
The vertical strip is for the world average 

r ee x r MM /r value. 



Figure [5] shows the comparison of our results with 
those of the previous experiments. The grey line 
shows PDG average and the error for the r ee xr MM /r 
product measurement. The new KEDR results are 
the most precise. Results are in good agreement 
with each other and with the world average value of 

Tee X / r. 
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Accounting the correlations in the r ee xr ee /r and 
r ee x r MM /r errors the mean value is 

r ee xr«/ T = 0.3320±0.0041 (stat.)±0.0050 (syst.) keV. 

With the assumption of lepton universality 
and using independent data on branching fraction 
B{J/ip -> e+e~) = (5.94 ±0.06)% leptonic and 



full widths of J/tp meson were determined: 

r« = 5.59±0.12keV 
r = 94.1 ± 2.7 keV 

These results are in good agreement with the world 
average [ifl and with the results from BaBar Q and 
CLEO-c [6] experiments. 
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